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Abstract 

Hybrid resonance is a physical mechanism for the heating of a magnetic 
plasma. In our context hybrid resonance is a solution of the time harmonic 
Maxwell's equations with smooth coefficients, where the dielectric tensor is 
a non diagonal hermitian matrix. The main part of this work is dedicated 
to the construction and analyzis of a mathematical solution of the hybrid 
resonance with the limit absorption principle. We prove that the limit 
solution is singular: it is constituted of a Dirac mass at the origin plus 
a principle value and a smooth square integrable function. The formula 
obtained for the plasma heating is directly related to the singularity. 

1 Introduction 

It is known in plasma physics that Maxwell's equation in the context of a strong 
background magnetic field may develop singular solutions even for smooth co- 
efficients. This is related to what is called the hybrid resonance [TTJ [TB], [7] 
for which we know no mathematical analysis. Hybrid resonance shows up in 
reflectometry experiments [141 113] and heating devices in fusion plasma [15] . 
The energy deposit may exceed by far the energy exchange which occurs in 
Landau damping ,16, 22] • In the physical literature it is somewhat considered 
as surprising: for example the 2x2 upper block of the dielectric tensor (|1.2p 
generates what is called X-mode equations, that is extraordinary equations, 
as opposed to the classical O-mode equation (O for Ordinary) related to the 
right-bottom coefficient in (|1.2[) . Hybrid resonance is a non damping dissipativc 
phenomenon: a singularity of the solution makes it stronger in some sense than 
the Landau damping; for example Landau damping needs kinetic effects to be 
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modeled, while hybrid resonance can be derived from the coupling of a fluid 
model with the non electrostatic part of the Maxwell's equation (derivation of 
the basic model is proposed in the appendix). Since the mathematical solution 
is not square integrable, hybrid resonance is also a paradoxal and non standard 
phenomenon in the context of the mathematical theory of Maxwell's equations 
for which we refer to [T2J [TUJ [3TJ [21] . The situation can be compared with the 
mathematial theory of metametarials. In [32] [33] the electric permittivity and 
magnetic susceptibility tensors are degenerate -i.e. they have zero eigenvalues- 
in surfaces, but they remain positive definite. In this case, the solutions are 
singular, but the problem remains coercive. In [5] [5] the coefficient changes in a 
discontinuous way from being positive to negative. In this situation coerciveness 
is lost, but as the absolute value of the coefficient is bounded below by a positive 
constant, the solutions are regular. In our case we have both difficulties at the 
same time. As the coefficient a (see below) goes from being positive to negative 
in a continuous way, its absolute value is zero at a point, and, in consequence, 
our problem is not coercive and there are singular solutions. 

Our purpose in work is to construct and analyze in Theorem 11.11 a mathe- 
matical solution with a hybrid resonance in slab (planar) geometry. We will use 
the limit absorption principle to construct the relevant solutions. We will show 
that an original singular integral equation is attached to the Fourier solution. 
Introduced in the seminal work of Hilbert []j|] and Picard [M] , this type of inte- 
gral equation is referred to as integral equation of the third kind, by comparison 
with the more classical equations of the first and second kind. Some references 
about this type of equations may be found in (3] [27] for mathematical analysis, 
and (29] El E] for relation with theory of particles or plasma physics. Our re- 
sults are reminiscent of those of Bart and Warnock [3] , even if our kernel does 
not satisfy exactly their hypothesis since it is less regular: that is the solution 
if the sum of a Dirac mass plus a principal value (plus a regular part). In their 
work it is stressed that non uniqueness is the rule for such equations. In our 
case, we are able to obtain uniqueness by means of the limit absortion principle 
which is a physically based selection principle. One originality of this work is 
the analysis of the properties of this singular equation for which we found no 
equivalent in the classical literature Q] E] . 

The model problem is based on the time-harmonic Maxwell's equations 



where the unknown E is the electric field, uj is the frequency and c the velocity 
of light. Plasmas with strong background magnetic fields can be represented by 
the cold plasma model [HI [11] for which the dielectric tensor is defined by 




(1.1) 



e{x) = P{x) 



(l 
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P{xf. 



(1.2) 



2 



A derivation of this dielectric tensor is given for convenience in our appendix. 
The incident wave with frequency L) is typically generated by an antenna on 
the boundary of the plasma. Such devices are actually being studied for the 
purposes of rcficctometry and heating of magnetic fusion plasmas in the context 
of the international ITER project: the ITER project is about the design of 
new Tokamak with enhanced fusion capabilities [3D]. The parameters of the 
dielectric tensor are the cyclotron frequency 

eB 

oj c = 

m e 

which is taken as constant in our work, the plasma frequency 



Som e 



which depends on the electronic density N e . The modulus of the background 
magnetic field Bq > will be assumed constant for simplicity in our work, e > 
is the absolute value of the charge of electron, m e is the mass of electrons and 
£o is the permittivity of vacuum. In this work we consider that the orientation 
of the background magnetic field is fixed, namely P is the identity matrix and is 
discarded in the equations. Boundary conditions can be of usual types, that is 
metallic condition n A E = 0, non homogeneous absorbing boundary condition 
like curlE+iXnAE = g on some parts of the boundary or even natural absorbing 
boundary condition at infinity. 

The equations for the transverse electric (TE) mode, E = (E x ,E y ,0), and 
E x ,E y , independent of z, which correspond to the 2x2 upper- left block in (|1.2j) 
are 

W +d y E x ~d x E y =0, 
d y W -a^E x -iS^Ey = 0, (1.3) 
-d x W +i5<E x -a<E v =0, 



where the coefficients are 
a = 1 



^1 oj c u>1 



w(w 2 -ui 2 c )' 



and W is the vorticity. 

In the plasma community this system is referred to as the X-mode equations, 
where the letter X stands for extraordinary mode or extraordinary waves. We 
suspect the reason is the non standard behavior of the solutions of this system. 
The case where to = uj c , i.e., when the frequency of the incident wave, w, is 
equal to the cyclotron frequency, uj c , will not be considered in this work. That 
is we consider that oj 7^ w c . If w < lo c it is called a low hybrid resonance. The 
other case u> > oj c is denoted as the upper hybrid resonance. On the other hand 
we will assume that the diagonal coefficient a is smooth and vanishes at x = 0. 
This configuration corresponds to the hybrid resonance. 
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x=-L 

Figure 1: X-mode in slab geometry: the domain. In a real physical device an 
antenna is on the wall on the left and send an incident electromagnetic wave 
through a medium which is assumed infinite for simplicity. The incident wave 
generates a reflected wave. We will characterize the antenna by the knowledge 
of the non homogeneous boundary condition (|I.4j) . The medium is filled with a 
plasma with dielectric tensor given by (|I.2j) . 

To be more specific we consider the simplified 2D domain 
= {{x,y) e M 2 , —L < x, yeM., L > 0} . 

We supplement the X-modc equations (|I.3[) with a non homogeneous boundary 
condition 

W + i\n x E y = g on the left boundary x = —L, A > 0, (1.4) 

which models a given source, typically an radiating antenna. In real Tokamaks 
this antenna is used to heat or to probe the plasma. We will consider slab 
geometry also for the sake of simplicity, that is all coefficients a and 8 arc 
functions only of the variable x 

d y a = d y 5 = 0. 

Other assumptions which correspond to the physical context of idealized rcflcc- 
tomctry or heating devices are the following. We assume that 

S eC^-L, oo[, < 6- < 5(x) < 8+, Vx. (1.5) 

The coefficient a vanishes at a finite number of points and a(x) = implies that 
a'(x) 0. For the sake of simplicity we will now suppose that a only vanishes 
at x = 0. More precisely 

a e C 2 [-L, oo[, a(0) = 0, a'(0) < 0, (HI) 
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Figure 2: X-mode equations in slab geometry: the physical parameters. The 
electronic density x <— > N e (x) is low at the boundary, and increases towards a 
plateau. The background magnetic field Bq is taken as constant for simplicity. 



and 

a_ < a(x) < a+, Va; e [— L,oo[, and < r < 

(H2) 

where H > 0. We will also assume that the coefficients are constant at large 
scale: there exists 5^ and so that 

5{x) = Sqo and a{x) = H < x < oo. (H3) 

We assume that the problem is coercive at infinity 

a 2 0O -S 2 O0 > 0. (H4) 

An additional condition is defined by 

4\\S\\loH < r. (H5) 

It expresses the fact that the length of the transition zone between x = and 
x = H is small with respect to the other parameters of the problem. One 
can refer to figures [5] and Q] for a graphical representation. This hypothesis is 
physically very reasonnable. Assuming (|H1|I - (|H4[) . it is known in the physical 
community that the problem is highly singular at the origin. Our main result 
can be summarized as follows. We denote by g the Fourier transform of g, 

9(0) ■= [ g(y)e- W ydy. 

We need the uniform transvcrsality assumption IH6I which is a generalization of 
assumption IH5I See Section [7J 

Theorem 1.1. Assuming HHliH6\) , for every g € L 2 (M.) with g of compact 
support there exists a solution in the sense of distributions of lil.2\) with boundary 



^ , \fxe[-L,H] 
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Figure 3: X-mode equations in slab geometry: parameters of the dielectric 
tensor deduced from the value of the physical parameters described in figure [5J 
assuming that lo > uj c . The coefficient a decreases from positives to negative 
values. It crosses the axis with a slope bounded from below by r. The coefficient 
8 is positive and bounded. Since the electromagnetic wave is strongly absorbed 
for x > H, we simplify by taking all coefficients constant for x > H because it 
does not change the physics of the problem. 



condition \1.J$ that goes to zero at infinity. Moreover, unless the source term 
g is identically zero, the electric field E x does not belong to ^J oc ((~ L,oo) x 
(—00,00)). The other components E y and W are more regular, they belong to 
L 2 {{— L,oo) x (—00,00)). 

Remark 1. The exact form of the solution is provided in |7.5[ j, where we have 
taken uj = c = 1 . The energy lost by the electromagnetic wave, which is also the 
one transmitted to the plasma, is given in formula |7.6| ) which is the Fourier 
characterization of related physical formulas provided at the end of the appendix. 

The result will be obtained with the limit absorption principle combined with 
a specific original integral representation of the solution. The loss of regularity 
of the electric field is counter intuitive with respect to the standard theory of 
existence and uniqueness for solutions of time harmonic Maxwell's equations 
[T2l [TOl [2T1 l3Tj . The essential part of the proof consists in showing that the 
Fourier transform E x may be composed of three contributions: a Dirac mass at 
x = 0; a non integrable function proportional to ^jy, that is interpreted as a 
distribution in the sense of principal value; and a regular part. The condition 
o/(0) < guarantees that the coefficient in front of the Dirac mass is finite. 
Moreover, the condition (|H5[) simplifies some parts of the mathematical analysis. 
The solution is a priori non unique since the limit absorption principle generates 
two solutions depending on the sign of the regularization. The heating of the 
plasma (|7.6j) is directly related to the singular part of the solution. 

This work is organized as follows. Section [5] is devoted to basic considcr- 
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ations. In the next section we introduce a regularization parameter, and we 
propose a specific integral representation of the solution. After that we recall 
the Plemelj-Privalov theorem and explain why it cannot be used directly for 
our problem. Section [5] is where we prove the properties of the solutions of 
the regularized equations. In particular, we show that one basis function has 
a fundamental singularity. Next, in Section [5] we define the limit spaces. The 
main theorem is finally proved in section [7] Additional material is provided in 
the appendix. 

2 Basic considerations 

In this section we rederive the phase velocity, compute the analytic solutions of 
the simplified Budden problem and introduce the limit absorption principle. 

2.1 Phase velocity- 
Recall that the phase velocity measures the velocity of individual Fourier modes. 



2.1.1 Constant coefficients 

Let us consider first that a and S are constant at least locally. A plane wave 
(E x , E y ) = Re^ x+k ^-y\ R G C 2 , is solution of X-mode equations (JOJ) if and 
only if 



k\ 
-hi k 2 



-fci fc 2 
kl 



a i6 
—i5 a 



R = 0, 



k = (ki, k 2 ) G 



We assume that c = 1 for simplicity. We set k = \k\d with d = (cos 9, sin 9) the 
direction of the wave. The phase velocity v v = ^ is solution of the eigenvalue 
problem 

sin# 2 — v^a — cos 9 sin t 
- cos 9 sin 9 + iv 2 „S 



cos# 2 — v^a 



R = 0. 



The determinant of the matrix is 



D = v% (a 2 -5 2 



i^- 



Sctting D = we obtain the phase velocity, 

9 a 



5 2 ' 



2.1.2 Non constant coefficients 

Let us assume for example that a = —x and that 5 — 1 which is locally compati- 
ble with the general assumptions. This represented in figureSJ We plot in figure 
[Hthe phase velocity as a function of the horizontal space coordinate. When the 
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phase velocity is real we are in a propagating region, and when the phase ve- 
locity is pure imaginary we are in a non-propagating region. One distinguishes 
two cutoffs where the local phase velocity is infinite 

Cutoff : a(x) = ±S(x) 

and one resonance where the phase velocity is null 

Resonance : a(x) = 0. 

This structure is characteristic of the hybrid resonance. 



V 



propagativc 
vl > « 


tran; 

3 V 2 = 

<f> 


rition 
= o 
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non propagativc 




x = 


= -1 X 


= : 


: = 1 



Figure 4: Sign of the square of the phase velocity. In this example a = —x and 
S = 1, so that v 2 = j^j. 

Remark 2. In what follows we always take ui = c = 1. 
2.2 The Budden problem 

In the case where the solution is independent of y, what for the plane waves 
corresponds to normal incidence, that is 6 = 0, the system (| 1 . 3[) is called the 
Budden problem [11] 

W - E' y =0, 
—aE x — iSEy = 0, 

-W + iSE x -aE y =0. 

After elimination of E x and W we obtain that, 

This equation can be solved analytically in some cases which helps a lot to 
understand the singularity of the general problem. Let us consider that a = — x 
and 5 is solution of — — x = — j + j:- The positive solution is 



S(x) = x jx 2 - | + 1 > 0. 
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The y-component of the electric field is solution of 

1 1 



E y + [-l + x) E y = - (2 - 1} 

This equation is of Whittaker type [2 [2] ■ It is a particular case of the confluent 
hypergcometric equation, and can also be rewritten under the Kummcr form. 
The 1 general theory shows that the first fundamental solution is regular 

v(x) = e~^x 

Indeed v'(x) = e~i (l — |) and v"(x) = ( — 1 + f), so that 

-- + - )v = 0. 
4 x J 

Let us consider a second solution w with linear independence with respect to 
the first one. The linear independence can be characterized by the normalized 
Wronskian relation 

v(x)w'(x) — v'(x)w(x) = 1. 
Seeking for a representation w = vz, one gets that 

/dx f dx i' f c"^ 

— => w = v / — = xe~ x/2 / — . 
v 2 J ir J x l 

Moreover, HE], 

% = --+/ - = --+E i (x) 

X^ X J X X 

where Ei{x) is the Exponential-integral function. It follows that, 

u(x) = -e x ^+xe- x/2 Ei(x). 

Furthermore [18] , 

OO j 

Ei(x)=]n\x\ + \Y, Ay- 

3=1 3 J ' 

It follows that, 

u{x)=-l + xh\\x\+0{\x\), \x\^Q. 

We notice that the second function w is bounded, but non regular at origin. 
It shows the subtleties associated with the singular Whittaker equation (12. ip . 
Nevertheless we note that the general form of the y component of the electric 
field of the Budden problem is bounded 

E y = av + bw => E y e L°°] - e, e[. 
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The x component of the electric field is more singular. It is a linear combination 
of two functions, the first one which is regular and bounded 



El(x) = ,V -i «(*) = ie-*^a — - + 1, 

and the second one singular which is singular at origin since w(0) = 1 



The general form of the x component of the electric field is a linear combination 
of these two functions. Since E™ £ L 2 ] — e,e[, we notice that the electric field 
is not a square integrable function in general. 



2.3 Limit absorption principle 

We will develop a regularized approach to give a rigorous meaning to the solution 
at all incidences. This regularized approach is based on the limit absorption 
principle. One considers a parameter \i =/= (the precise sign will be justified 
later) and the regularized problem with unknown (Eg, Ely, W 11 ) 

W» +d y E% -d x Eg = 0, 

d y W^ -(a{x)+i/j,)E% -iS(x)E^ = 0, (2.2) 
-d x W» +iS(x)E% -(a(x) +in)E y l =0. 

The regularization parameter \i can be interpreted as a small collision frequency. 
See Appendix B. 

A further simplification consists in Fourier reduction. Since the coefficients 
do not depend on the y variable, one can perform the usual one dimension 
reduction. The system that will be studied in this article is obtained by applying 
the Fourier transformto the regularized system (|2.2[) . Denoting the unknowns 
(U, V, W) it yields 

W +i6U -V =0, 

i6W -(a(x) +ifj,)U -iS(x)V =0, (2.3) 
—W +i5(x)U -(a(x)+ifj,)V =0. 

Here the notation ' denotes the derivative with respect to the x variable. 



3 A general integral representation 

We begin by some notations. Let us denote by (A^,^) the two fundamental 
solutions of the modified equation 

- u" - (a(x) + ifj,)u = 0, (3.1) 
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with the usual normalization 

A„(0) = 1, 4(0) = andS M (0)=0, = 1. (3.2) 

Various continuity estimates of and i? M are derived in the appendix for the 
sake of the completeness of this work. Let us denote T> 6 Z the operator i0d z — i8(z) 
applied to any function h, that is 

V e z h = i6d z h - iS(z)h. (3.3) 

Let us define the kernel 

k"(x, z) = B M (z)A M (ar) - B^A^z). (3.4) 

Next we define 

!V e V k ti (x z) 
x /\ [ ' \ for G < z < x or x < z < G, 

0, in all other cases. 

Let us define the kernel sequence by 



(3.5) 



K%(x,z;G) = U * V **} X : t) K^(t,z;G)dt. (3.6) 
Jg a ( x ) + W 

The sum is 

oo 

K e >"{x, z;G) = J2 K^(x, z; G). (3.7) 
n=0 

The integration domain is centered on G, that is 

supp (ifi'^O, •; GO) C {(x,y) e R 2 ; G < z < x or a; < z < G} = V G , (3.8) 
which yields as well 

supp (#*•"(•, -;G)) c% 

Proposition 3.1. Any triplet (U, V, W) solution of the regularized system (|2.3[) 
admits the following integral representation. 

• One first chooses an arbitrary reference point G £ [— L, oo[. 

• The x component of the electric field is solution of the integral equation 

/• BWVC,,,) f^. (3.9, 

where the right hand side is 

F e ^(x) = a G V 9 x A^x) + bcVlB^x) (3.10) 
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and the kernel is given in \3. $\3.J$ . The solution of this integral equation 
is naturally provided by the resolvent integral formula 



U(x) = . . ~\~ I K,' i '> 1 (x,z;G) , b . v ' dz (3.11) 



where the resolvent kernel is constructed in (3^ 
• The y component of the electric field is recovered as 



V{x) = acA^x) +b G B fl (x) + / V° z k^{x,z)U{z)dz, (3.12) 
and the vorticity is recovered as 

W(x) = a G A'^{x) + bcB'^x) + I d x V s z k»(x, z)U(z)dz. (3.13) 



The two complex numbers (og,6g) solve the linear system 

a G A li (G) + b G B^(G) = V(G), 
a G A' fl (G) + b G B' )1 (G) = W(G). ^' i4j 

Proof. Eliminating W from the first and third equations of (|2.3[) gives 

-V" - (a + ifj)V = f with / = -iOU' - iSU. 

Since the Wronskian is constant, it follows from the normalization (|3.2p that 
A^B'^ — A'^B^ = 1. Then, from the variation of constants formula, 

V(x) = ajA^x) + bfB^x) + [ f(z)k' i (x, z)dz, Vx. (3.15) 

JG 

where a/ and 6/ are two integration constants. Now we replace / by the corre- 
sponding function of U and perform the integration by part 

U'(z)k fl (x,z)dz = U{x)k fl {x,x) - U(G)k fl (x,G) - / U{z)d z k^{x,z)dz. 



Since k^{x,x) = there is a simplification. Therefore (|3. 15|) yields (|3.12[) with 
a G = a f + i6U(G)Bf,{G) and b G = b f - i6U{G)A^(G). Next we eliminate W 
from the first and second equations of (|2.3p and obtain 



-i6V' -6 2 U +(a + in)U + i5V = 0. (3.16) 
The derivative of (|3 . 12[) yields 

V'{x) = a G A' ll {x) + b G B^{x)+ f d x V e z k tl {x 1 z)U{z)dz + V e z k"(x,x)U{x). 
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Since T> 9 z k' x {x,x) = id (A^B'^ - B^A'^) = i9, one gets the identity 

V\x) = a G A'(x) + b G B'^{x) + / d x V e z k^( Xl z)U{z)dz + i6U{x). 

Jg 

Plugging this expression in (|3.16[) and performing all simplihcations we obtain 
the integral equation (|3.9[) . Finally, we get the last integral formula (|3.13p 
from W = —i8U + V . The linear system p,14[) is obvious from (|3. 12113.1^) at 
x = G. □ 

Following [53] , the equation (|3.9p is an integral equation of the third kind in 
the case fi = 0. In this case the theory is rather incomplete regarding existence 
and uniqueness [3]. However as long as fi ^ 0, the solution based on these 
integral equations is uniquely defined. Then, the question is to determine the 
behavior of these solutions when fi goes to 0. Moreover, different choices of G 
will give different kind of information. A strategy to study of the limit solution 
— > can be the following: Choose an optimal G, so that a) the integration 
constants (ac, &g) a^e easy to determine, and b) the resolvent kernel K, •; G) 
admits a limit as fi — > 0. Considering the form of the right hand side in p. lip , a 
convenient tool is the Plemelj-Privalov Theorem (231 [26] . Unfortunately, we will 
see that a fundamental singularity of the kernel /C e,A1 (-, ■; G) prevents any simple 
limit procedure. A more convenient technique will be proposed in Section [5] 



4 Singularity of the kernels 

A fundamental tool in order to pass to the limit in singular integrals is the 
Plemelj-Privalov theorem [23l [26]. However, to apply this theorem to pass to 
the limit fi — > in equation p. lip it is necessary that the kernel /C e,AI (x, z) be a 
Holder continuous function of z for each fixed x. Unfortunately, this regularity 
is not available in our case. To illustrate this phenomenon, we study only the 
first term of the series (|3.6p that defines K. 6 '^, namely 

fc^(z,z):= a ~\ [ ' K (4.1) 

We consider two cases. 



4.1 First case: G ^ 

In this case there exists (0, z) £ T>q with z ^ 0. In the limit case /i = one has 
that /Ci°(x,z) admits the local expansion: 

K, e 1 °{x 1 z)^-^ ) V s x Vlk^x,z). 

Therefore, Kf{°{x, z) blows up as x — > 0. 
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4.2 Second case: G = 

We turn to the case G = 0. We begin with a preliminary result. 
Proposition 4.1. One has 

{Vlvlk^ix, x)=0 Viel. (4.2) 
Proof. Indeed by construction 

(V e x V e z k^)(x,x) = -5(x)S{x)k' i {x,x) 

+0 S (x) ( (d x ft" ) (x, x) + (d z k» ) {x, x) ) - e 2 (d x d z W ) {x, x) . 

We notice that by definition k^(x,x) — for all x so the first contribution 
vanishes in (T) x 'D z *k tl ){x,x). One also has that 

(d x k»)(x,x) + (d z k> 1 )(x,x) 

= B^A'^x) - B'^x)A^(x) + B^x)A^(x) - B^x)A'„{x) = 0, 
so, the second contribution vanishes also. Furthermore, 

(d x d z k»)(x,x) = B'^A'^x) - B^(x)A^x) = 0. 

This completes the proof of equation (|4.2j) . □ 

Proposition 4.2. The limit kernel 

V x V e z k^=°(x,z) 
a(x) 

belongs to L°° (D ). 

Proof. A first order Taylor expansion of r D x r D e z k tl around yields 

V e x V 9 z k»(x, z) = a^x + p^z + Q(\x\ 2 + \z\ 2 ). 

Notice that (|4.2|) implies /3 M = — a M . The coefficient a M is easily computed using 

(2>«Z>f *")(*, 0) = V x A,(x)V e z B,(Q) V e x B,(x)V e z A,(Q) 

and the definition (|3.UI3.2j) . One gets that V e x A^(x) = -iS(0) - iS'(0)x + 9/j.x + 
0{x 2 ) and V e x B^(x) = id - iS(0)x + 0{x 2 ). So 

(V 9 x V e z k")(x, 0) = (-id(0) - iS'(0)x + 0/juc + 0{x 2 )) id 

- (i6 - iS(0)x + 0(x 2 )) (-id(0)) 
= (<S(0) 2 + 66' (0) + i6 2 n) x + 0{x 2 ). 
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This coefficient a M being constant, one obtains that 



¥x{z) 



V e x V e z k^°(x, z) _ (S(0) 2 + 95' (0)) (x-z) + Q{\x\ 2 + \z\ 2 ) 



a(x) 



a(x) 



(4.3) 



This expansion is valid for (x, z) G T>o (the domain T>o is defined in (|3.8|) 1: in 
this case \x — z\ < \x\ and \z\ < \x\. Moreover, since a(x) = x(a'(0)x + 0{x)) 
we obtain that, 

\8(0) 2 + 65'(0)\ 



1^(^)1 < 



0(\x\). 



Since there is no such difficulty for x away from 0, this inequality ends the proof 
of the proposition. □ 

Remark 3. A similar property holds for T>x ^^)+i^ w hich also belongs to 
L°° (T>q) for all 9 and uniformly for fi G [—1, 1] \ {0} ; that is 



V e x V e z W(x,z) 



a{x) + ifi 



<C e , MG [-1,1]\{0}. 



(4.4) 



Such estimate is sufficient to control some L°° bounds of the series that 
defines the iterated kernel K, e,fJ, (x, z; 0). However, L°° bounds are not sufficient 
to show that K, e,li {x, z; 0) is of Holder class in z in the vicinity of x = : That 
is, one cannot pass to the limit using the Plemelj-Privalov theorem for all values 
of the parameters involved in (|3.61 13. 1 1[) . This is why we will develop another 
approach to give a meaning to the limit value. 



5 The space (/i ^ 0) 

The solutions of the integral equations evidently belong to a vectorial space 
of dimension two: see also (|5.2[) . In a first stage we will design a particular 
basis in this space, in a second stage we will study the properties of the two 
basis functions. A careful analysis of this singularity will allow to show that 
one basis function (more precisely the x component of electric field) is the sum 
of a singular part a ^ +i ^ plus a term which is bounded in LP (1 < p < oo) 
uniformly with respect to /i. It will be the central result of this part. 

For the simplicity of notations, we restrict the parameter to < n < 1 with- 
out loss of generality. The extension to negative [i will be considered in section 
(|6.2[) . We define the vectorial space of all solutions of the X-mode equations 

X 6 '^ = {x i ^ (U(x), V{x),W{x)) , for all solutions of the system (gS])} . 

(5.1) 

One may also use the notation 

This section is devoted to the analysis of this space. 
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Remark 4. The property that dim X e,M = 2 is also evident considering the right 
hand side of the integral equation HS.9\) . 



By elimination JJ ^ in (|2.3|) . one gets a system of two coupled ordinary 
differential equations 

I { ) - **w ( v w z ) m 



with 



ggfg) 



A W ~ 8(xf / \ 98 { x) l 

\ a(x)-\-ifi ^ J ' a(x)+ifi / 

In the case /i 7^ the matrix is non singular for all x, which gives a meaning to 
the regularized problem. One notices the matrix is singular for /j = 0. 

Lemma 5.1. Take two solutions (V 6 ^ ,W e '») and (y e <», W 6 '^ of (h*lfy. De- 
fine the Wronskian 

W(x) = V e ^(x)W e 'V(x) - W e ^{x)V 9 '^(x). (5.4) 

Then the Wronskian is constant 

W(x) = W(0), Vx. 

Proof. The system (|5.2[) main be rewritten as 



_d ( V \ ( a b \ f V 
dx V W J ~ \ c -a )\W 

Therefore 

-T-W = 4- (V^W^x) - W^V^x)) 
dx dx \ I 

= (aV + bW) W + v(cV- aW) - (cV - aW) V - W (aV + bWj = 
since all terms cancel each other. □ 

5.1 The first basis function 

Next we desire to particularize a convenient basis in this space. The first basis 
function 

= (u^, W^) e (5.5) 



is the natural one which is smooth at the origin: U 1 ,M (0) = 0. For that reason G 
is chosen to be the origin in this subsection, so that the corresponding integral 
equation has a bounded right-hand side and a bounded kernel. It is naturally 
characterized by 

^/' M (0) = i6», and wf' M (0) = iS(0) (^ 0). (5.6) 
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Proposition 5.1. The basis function i5.5\) is uniformly bounded with respect 
to \i: for any interval 9 £ [#_,#+] and any H €.]L, oo[, there exists a constant 
independent of fi such that 





+ 


ye* 




+ 






L°^(-L,H) 






(-L.H) 





L°^{-LM) 

Proof. The right hand side in the integral equation (|3.9[) is 



< C. 



(5.7) 



with h^(x) = iOVlA^x) + i6(0)V%Bp(x). 



a(x) + ifi 
With the choice (|5.5I) one has 

h»(O)=iO(-i6(O))+i6(O)(i0) = O V/i. 
Therefore the right hand side of the integral equation 

a[x) + ifj, 

is bounded in L°°(—L, H)) uniformly with respect to [i. The solution [/f (|3.1ip 
is also bounded, since by the results of Subsection 4.2 the kernel K '^(x, z, 0) 
is also uniformly bounded. These bounds are uniform with respect to /i. The 
integral representation (|3.12p of the V 1 ' M yields that V x ' M is also bounded. It 
is similar concerning the integral representation p. 131) of the Wf' M , so W 1 ,M is 
also bounded. □ 



5.2 Behavior at infinity 



Hypothesis (|H3[) allows to study a simplified model with constant coefficients for 
x > H . In fact, it corresponds to a system as in (|5.2[) with constant coefficients, 
which matrix will be denoted A 6 ^. 

Proposition 5.2. The matrix A^t 1 has two distinct eigenvalues. The first 
eigenvalue \ e4i has a positive real part. The second eigenvalue is — A 8 '^. 

Proof. The eigenvalues are solution to the characteristic equation 
A 2 -tr(A^)A + dct(A^)=0 

where tr(A^) = and 

51 



det(A^) = «oo + in 
The real part is 

real(det(A^))=a oo -0 2 -4^ 



S 2 

OO 



"So + V \ a lo + I- 1 - 2 

and is therefore negative due to the coercivity assumption (|H4|) . So the usual 



square root \ 04i = y — dct(A 00 AI ) has a positive real part. The other one has a 
negative real part. □ 
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As a consequence any U £ X"' 11 is at large scale a linear combination of the 
exponential increasing function and a exponential decreasing function 

U(x) =c + R + e x °'' Lx + c-R-e- xe '* x H<x (5.8) 

where R + € C 3 and R_ G C 3 are constant vectors and (c + , c_) G C 2 arc 
arbitrary complex numbers. Regarding the structure of the matrix and using 
the second equation of the system (|2.3|) . one gets that R + = (r+, r\, r\) with 

x i6r\ iS{H)r\ 2 _ 3 _ / TT, ^ d5{H) 



The other vector i?_ = (ri. , r 2 , r 3 ) is characterized by 



rt = 1 777^— - ri = -V-det(^) 



a(H) + ifi a(H)+ifi * a(H)+ifj, 

One notices that i?+ and i?_ are well defined for all /i £ K, in particular even 
for /X = 0. 



Proposition 5.3. The first basis function \5.5\) is exponentially growing at large 
scale (/j, 0). 

Proof. For the sake of simplicity, denote U 1 '' i = (Ux, Vx, Wx), dropping the 8s 
and fis. Then from system (|2.3p one gets 

Wi +i8Ux —V{ = 0, 

idWx -(a + ifx)Ui -iSVi = 0, 
-W{ +i5Ui -(a + ifi)Vi = 0. 

Multiplying the second equation by Ux and the third one by V\ , the sum writes 

iOWjh-WlVi- (a\Ui\ 2 + a\Vi\ 2 + iSVxTh - i8UxVx)-in (\Ux\ 2 + \Vi\ 2 ) = 0. 

On the other hand an integration in the interval }M, N[ yields 



N 
M 



/ (iOWjh - W{W) dx 

JM 

iQWjJl + WxV?} dx - W x {N)Vi{N) + Wi(M)Vx(M) 

N 

Wi\ 2 dx - Wx{N)Vx(N) + Wx(M)Vx(M), 



<M 

where we used the first equation. We obtain the identity, 

/ (\Wx\ 2 - a\Ux\ 2 - a\Vx\ 2 -iSVxVx+iSUxVi) dx~ifi / (\Ux\ 2 + \Vx\ 2 ) dx 
Jm Jm 

(5.9) 
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= W x (N)Vi(N) - Wi(M)Vi(Af). 
Splitting between the real and imaginary parts, one gets the important relation 

/■N 

M/ (l^i I 2 + \Vi I 2 ) dx = Im (Wi (M)Vi(M)) -Jm(Wi(N)Vi(N)) (5.10) 

J M 

which is true in fact for any element in X e,M and for any M < N. 

Let us take M = 0: so Vi(0) = y^yWi(O) and Im (Wi {M)V[ (M)) = 0. 
Therefore 

/•iV 

/* / (iC/il 2 + l^l 2 ) dx = -Im(Wi(JV)T^(JV)) . 

It shows that Wi(iV)Vi(-/V) for iV — > oo. In other words the first basis 
function does not decrease exponentially at infinity. Considering (|5.8[) it means 
that this function is exponentially increasing at infinity. □ 

5.3 The second basis function 

The second basis function 

is built with two requirements. 

• It is exponentially decreasing at infinity, that is 

U^(x) =c_i?_e- Afl '^, H<x, (5.11) 

for some c_ G C. 

• Its value at the origin is normalized with the requirement 

ifxU^iO) = 1. (5.12) 

To ensure that these conditions can be satisfied, consider the third function 

U|'" = (^l K?* W^)(x) = R^e- xa " x H<x, (5.13) 

where R- and A_ are defined in section I5~2l smoothly extended so that Ug'^ 1 € 
X e 'f". The identity 

n J m (l^'T + l^l 2 )** 

= Im (V 3 e ' M (A/)"t^p(A/)) - Im (w^ \N)vf~» '(N)} 
with N — > oo and M = shows that 

A 4 ^ (Va'f + \V°'r) dx = Im "(0)^(0)) . 
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However, from if '"(O) = ^jW 3 (0) - j^jE/^O), so one gets 

(\U^\ 2 + |lf >f ) = - J^Im (wf ^(0)[f^(0) 



It shows that /7g' M (0) ^ 0. This is why it is always possible to renormalize with 
a parameter 

= C - U 3" C = ( 5 - 14 ) 

so as to enforce (|5.12p . 

Proposition 5.4. Wi£/i £/ie normalizations \5. 6)) and Il5.11]\5.l!fy . the Wron- 
skian relation takes the form 

V^{x)W^{x)-W^{x)V^{x) = 1 Vx. (5.15) 

Proof. It is sufficient to compute it at the origin 

if '"(0)Wf '"(0) - wf '"(0)lf '"(0) = iflWf '"(0) - ^(0)lf •"(O) 

= (a(0) + «/i)C/ 2 e ' M (0) = ifxU^(0) = 1 
using (|23]1 and thanks to (|5~T2|) . □ 

Remark 5. TTie va^we of the Wronskian L5.15]) is independent of \i. It will be 
of major interest in the limit regime fi — > 0. 

The non zero Wronskian shows (|5. 15|) shows that the two basis function are 
linearly independent. So they span the whole space 



Spanju^U^ 1 }, j u>0. 



5.4 Passing to the limit \i — > 

We now study the limit \i — > 0. An important result is that the first basis 
function admits a limit which is defined as a continuous function in C°[—L, oo[ 
and is independent of the sign of \i. On the other hand the second basis function 
admits a limit which is singular at x = 0. Moreover the limit is different for 
(i — > + and for fi — > 0~. The linear independence of these limits will be 
establish with a transvcrsality condition. 

5.4.1 The first basis function 

There is no difficulty for this case which is easily treated passing to the limit in 
the integral equation p. lip , choosing G = 0. The limit basis function is referred 
to as 

U? = (Etf,lf,Wf) 

is and will be called the regular solution by analogy with the terminology 
in scattering on the half-line. 
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Proposition 5.5. The first basis function satisfies for fj, — > 



-> 0. 

(L~(-L,_f/)) 3 



Proof. Setting G = both the resolvent kernel and the function F , . f 3 ^ are 
bounded for all /Lt 7^ 0, and admit a continuous limit with respect to fj, — > 0. The 
limit basis function is defined by 

Vf (0) = i0, and Wf(0) = iS(0) 0). (5.16) 

Then C/f is the solution of the non singular integral equation 



vim - r M*)*z - i ^m±mp^Mi j 

Jo a K x ) a \ x ) 

with k(x,z) = A(x)B(z) — B(x)A(z). This equation is well posed since the 
kernel is locally bounded and the right-hand side is a locally bounded function. 
The other components are 

V?(x) = i6A{x) + id{0)B(x) + f V 9 z k(x, z)U 9 l {z)dz 

Jo 

and 

W?(x) = iQA'{x) + iS(0)B'(x) + [ d x V 9 z k{x, z)\j{{z)dz. 

Jo 

Therefore the convergence is granted by classical dominated convergence theo- 
rems. □ 

The next result establishes that is still exponentially increasing at infinity 
with a technical condition. 

Proposition 5.6. Assume hypothesis (|H5[) . Then\J 9=0 increases exponentially 
at infinity. 

Remark 6. The constant 4 in the condition (|H5[) is probably non optimal. 

Proof. We drop the super-index • to simplify: that is (Ui, Vi, W\) stands for 
(Ui, V] , Wi). Let us consider the identity (|5.9j) which holds true at the limit 
[X = 

N 

(\W!\ 2 - a\Ui\ 2 - a\V x \ 2 ~ iSViTh + iSUiW) dx 

= W 1 {N)V[(N) -Wi(0)V[(0), 0<N<od. 

Since we consider the case 9 = 0, Vi(0) = 0. Notice also that W\ = V[, so the 
relation is rewritten as 

/ {\V{\ 2 - a\Ux\ 2 - a\Vi\ 2 - iSViUi + iSUiVZ) dx = W^N^N). 
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Let us proceed by contradiction: we assume that the function is exponentially 
decreasing at infinity. It yields 

/•OO 

/ (\V{\ 2 - a|J7i| 2 - a\V!\ 2 - iSVjh + iSUi W) dx = 0. 
Jo 

Notice that -a|f7i| 2 - a|F| 2 - iSVjh + iSUiVi > for x > H due to the 
coercivity property (|H4]I . Therefore it implies that 

r H — 

/ (IF/I 2 - ct|C/i| 2 - a|^i| 2 - WViC/i + ^C/iVi) < 0. 
Jo 

Next observe that U-i = —i—V\ , so that 

Since Vi(0) = and a(x) sa o/(0)x with o/(0) < (see hypothesis IH1 [). it is 
convenient to notice the proximity with the famous Hardy inequality that we 
recall, 

f H ^f- <4 f" u'(x) 2 , ueH 1 ^,!!), u(0) = 0, u ? 0. (5.17) 
Jo ^ Jo 



Since, thanks to hypothesis ()H2|) . 



it yields the inequality 



l<5|l °° ^) iT |yi ' |2dE - 1 " ( |yi '' 2 1 ^ IT/ 12 ~' IT/ 12 



where we used (|H5|) . Therefore V\ vanishes on the interval [0, if"]. So JJ\ vanishes 
and W\ also vanishes on the interval which is not compatible with Wi(0) = 
iS{0) ^ 0. 

□ 

Proposition 5.7. There exists a maximal value ^thresh > such that: If 
H^Hoo if < 4r and \0\ < #thresh, then Uf increases exponentially at infinity. 

Let us denote by {U§(H),V£(H),W$(H)) the solution to (£3]) for x > 
that satisfies f|5 . 13[) with /i = 0. 

Proof. Let us consider the function 

,7(0) = - Wl(H)V 3 e (H) (5.18) 
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By definition 



(V S (H),W 3 (H)) = 1 - A JJSS. - J- aoo + e*+ 5 M e -V-dct^=»H. 



S' 2 

This vector is real and always non zero. Therefore the function 9 H> f(6) is well 
defined. This function naturally satisfies two properties 

• cr(0) ^ since (Vi,W±) is exponentially increasing by virtue of the 
previous property. Indeed er(0) = if and only if the functions x \-¥ 
(Vi(x), Wi(x)) and x H> (V^(x),W°(x)) are linearly dependent, which is 
not true. 

• the function a is continuous since the first basis function is continuous wit 
respect to 6. 

Therefore there exists an interval around in which a (6) is non zero, which 
in turn yields the fact that Uf is linearly independent of U3. Therefore XJf is 
exponentially increasing. □ 



5.4.2 The transversality condition 

Passing to the limit in the second basis function near the origin is involved. 
Indeed we expect that the limit [/f is such that £/| ~ — for some local constant 
C. Therefore the limit is singular and special care has to be provided to avoid 
any artifacts in the analysis. 

Let us define the special Wronskian between the first and third basis func- 
tions 

a{d,fj) = V?^{H)W^{H) - Wl'"(H)Vg' lt (H). 

It is the natural continuous extension with respect to \i of the function 6 1— > a(9). 
We rewrite (|5 . 14[) as 

Plugging this relation in the Wronskian (|5.15[) one gets that 

This function is continuous with respect to [i. Moreover the function defined 
in (|5.18|) satisfies a(6) = <r(9,0). The transversality condition is defined as the 
condition 

a{6) ^ 0. (5.19) 

If the transversality condition is not satisfied, that is <j(0) = 0, then by continuity 
|£ e ' M | — > 00 for fj, — > 0. If a(8) = 0, then the first basis function and the third 
function are linearly dependent at the limit /j, = 0. It is of course possible to 
develop the theory in this direction, but it seems to us less interesting. Therefore 
we will always assume the transversality condition^ from now on. We postpone 
to the appendix some aspects of the case where is not satisfied. 

lr The "transversality condition" is a sufficient condition of linear independence. 
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Proposition 5.8. Assume the transversality condition H5.19\ ). Then for all 
e > one has the limit 



— — if 



-> 0. 



(i°°[e,°o[) 3 



Proof. Evident. 



□ 



In order to show that the second basis function admits a continuous limit for 
x < 0, the strategy is to solve the integral equation Q3.9[) from G = H backward, 
and to show that fine estimates on the solution give knowledge of the limit even 
for x < 0. 



/ 



Condition at infinity 



x=0 



Figure 5: Schematic representation of the real part of the limit electric field 
of the second basis function U 2 , /i > 0. Here the transversality condition 
o~(9) ^ is satisfied, which turns into a singular behavior at the limit [i — > 0. 



5.4.3 Continuity estimates 

The integral equation p.9[) is singular at the limit. The whole problem comes 
form the singularity at x = 0. By comparison with the standard literature 
[2"51 12"31 13TH [T2l [2^1 [31 [57] we found no convenient mathematical tool to analyze 
its properties. That is why we develop in the following new continuity estimates 
with respect to the parameters of the problem. On this basis we will manage to 
pass to to the limit [i — > 0. 

Let us consider a general solution U = (U, V, W) € X e,A1 of the integral 
equation (|3.9[) with prescribed data in H under the form 

V(H) = a H and W(H) = b H - 
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Let us introduce the compact notation 

= \a H \ + \b H \. 

Our goal is to obtain some sharp continuity estimates on the solution U with 
respect to \\H\\. The main point is to bound the constants uniformly with 
respect to < fi < 1 which is hereafter taken positive for the simplicity of 
notation. The reference point can be different from H as well, but non equal 
to zero. Once these continuity estimates are proved, they will provide enough 
information to define the limit (i — > of the second basis function. 

Proposition 5.9. There exists a constant Cg with continuous dependence with 
respect to 6 such that 

\U(x)\ < _5jjL=||tf||, 0<x<H. (5.20) 
\Jr- L x i + /i z 

Proof. Let us consider 

70= sup ||A M || w i,«.( 0) jr) + sup \\Bh\\wi.°°(o,h) ) (II^IU + 

\0<fi<l 0<fi<l J 

The integral equation (|3.9[) with G = H implies that 



s/ r 2 x 2 + /i 2 J x y^r 2 x 2 + /i 2 

where we used (|H2|) . Since T>^.T> s z k{x,x) = for all x, there exists a constant /3g 
such that 



So 



and 



\V e x V e z k(x,z)\\ L ^ ]0 H[ < (3 e \x - z\ < p e z for < x < z. 
Vr 2 x 2 + fJ, 2 \U(x)\ < le \\H\\ + fig / z\U{z)\dz 

J x 



H 



rx\U{x)\ <~fg\\H\\+ (3 e z\U{z)\dz, < x < H. 



The Gronwall lemma is useful to study this inequality. Indeed let us set g(x) 
J H \zU{z)\dz, so that the previous inequality is rewritten as 

-rg'{x) < le \\H\\+(3gg(x). 

Therefore < 7e||-ff|| + rg'{x) + /3 e g(x), that is 

Q<le\\H\\e^ x +r(e^ x g{x))' . 
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Next we integrate on the interval [x, H] and use the fact that g(H) = by 
definition. It yields 



•-H 



0< 7fl ||F|| 



- re^ x g(x), 



that is 



Finally one checks that 



g(x)< 3 io\\H\\ 

pe 



(5.21) 



y/r 2 x*+v?\U(x)\ < -fe\\H\\ + I3 e g(x) < e^^-^ le \\H\\ 

which proves (|5.20|) . □ 

Next define ||0|| = |V(0)| + |W(0)|. 

Proposition 5.10. There exists a constant Cg with continuous dependence with 
respect to 9 such that 

||0|| < C e (l + |ln M |)||if||. (5.22) 

Proof. We adopt the same notations as above. The integral expression of V 
(I3.12P with G — H yields the inequality 



\V{G)\<le\\H\ 



V z k(O,z).U(O) 



dz 



We notice that V z k(O,z) = (V e z k(0,z) - V e z k(0, 0)) + V e z k(0, 0). Since 
V 6 z k(0, 0) = m x k(0, 0) - iSk(0, 0) = i9 

one gets 

\V 8 z k(0,z)-i6\ <w\z\ 
for some constant rje > 0. Which gives 



Wo)l<7ell#ll + W z\u(z)\dz+\e\ 



H 



H 



U(z)dz 



By (f5~20]) Q < Cg\\H\\, and moreover, 











R := 


[ U{z)dz 


< Cg\\H\\ 


I 




Jo 







-dz 



<C fl ||ff|||ln/x| 



This completes the proof for | V2 (0) | . The term | W(0)| is bounded with the same 
method starting from the integral (|3.13[) and using the identity d x V z 'k fl (x, x) = 
i5{z)A fl (z). □ 
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An interesting question is the following. Let us consider the integral equation 
(|3.9p with G = 0. That is the starting point of the integral is the singularity. 
One may wonder if a direct use of the Gronwall lemma may yield valuable 
estimates, or not. It appears that a pollution with log (i terms render the result 
of little interest. 

Consider firstly for simplicity < x. Then (|3.9|) with G = turns into 



\U(x)\<C e .J 1 " 11 +C ( X \U{z)\dz (5.23) 
^/r^x z + [L A Jo 

where we used (|4.4j) to bound the kernel. The constant Cg > is chosen large 
enough. Set h(x) = J Q \U(z)\dz so that 

h'(x)<C e - 11011 +C e h(x). 
y r z x z + [i A 

Since h(0) = the Gronwall lemma yields the inequality 



h(x) < C'g / , 1 ' , dz 

that is after integration (0 < x < H) \h(x)\ < Cg||0|| (l + |ln/i|), for some 
constant Cg > with continuous dependance with respect to 9. Considering 
the bound (|5.22|) and the symetry between < x and x < in the integral (|3.9[) 
(with G = 0) one obtains the estimate 



U(z)dz 



< C'g"\\H\\ (l + |ln/j|r, -L<x<H. (5.24) 



Going back to (|5.23|) which is easily generalized to x < 0, one gets 

\U(x)\ < C e ( 1 + 1 + | In m| J (1 + I In Ml) ll#H, -L < x < H. 
\ \Jr i x z + fi z ) 

(5.25) 

By comparison of (|5.20p and (|5.25[) , it is clear that this technique generates 
spurious terms of order log n for positive x. It spoils the possibility of having 
sharp estimates also for negative x. With this respect, the rest of this section 
is devoted to the obtention of various sharp inequalites which are free of such 
spurious terms. 
Let us define 

Q(U) = '"(H)W(H) - W^{H)V(H). (5.26) 

This quantity is the Wronskian of the current solution U against the first basis 
function. It is therefore independent of the position H which is used to evaluate 
Q(U). 
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Proposition 5.11. There exists a constant Cg with continuous dependence with 
respect to 9 and a continuous function fj, i— > with e(0) =0 such that 



L 2 (-L,H) 



ttQ(U) 2 
a'(0) 



< C e e(/x)||if| 



(5.27) 



Proof. We consider positive /i to simplify the notations. The proof is easily 
adapted for negative fi. 

Consider the integral equation (|3.9p with G = 0. One gets 



apP^4(g) + b V e x B(x) 
a(x) + ifi 



X -J-.0 



a(a;) + i/i 



U(z)dz. 



Here (ao, &o) are a priori different from (a//, bn)- Due to (|5.26p . the normaliza- 
tion of Ui and thanks to Lemma 15.11 one has that 

a V x A(0) + b o V x B(0) = Q(U). 



So the integral equation can be written as 

Q(U) 



U(x) 



a(x) + ifj, 

v ' 

Si 



(5.28) 



V x A(x)-V x A(0) V x B(x) - V x B(0) 
+ a , , , . h b - 



a(x) + ifi 



a{x) + ifi 



a(x) + i/j, 



U(z)dz. 



S2 S3 
• The L 2 norm of the first term Si depends upon the value of 



D„ 



/' 



-L Oi(x) 2 + n 2 



dx. 



Make the change of variable x = fiw so that = J_ f ' t b f^yi +1 dw and 



Using the hypothesis (|H2p one has that |& M (u>)| > rw, r > 



dw 



0. Since J M 

minai 

Considering that 



^ < 00 and the point-wise limit of b^(x) is a'(0), the 
Lebesgue dominated convergence theorem states that lim + = rpTgyi ■ 



\Q(V)\<Cl\\H\\ 



(5.29) 



using (|5.26j) . there exists a continuous function (j, n- e 1 ^) with e 1 (0) = 
such that 



MllSil 12 



L 2 (-L-,H) 



7TQ(U) 5 



a'(0) 



< cy^)\\H\\ 



(5.30) 
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. The functions g^Mzg^g) and gi^MzgiM can bounded in L oo 
uniformly with respect to /z. So 



J H L \S 2 (z)\ 2 dz<c 2 e \ 



Estimate ([Q2"]) yields 



mII^II^^.^^qVi + IIumI) 2 !!^!! 2 



for some constant Cg > 0. 
• The last term 1S3 is 



\SM\ 



a(a;) + 



since the kernel is bounded (|4.4[) with respect to 6* and uniformly for 
fi G [0,1]. Inequality flQ3) implies that |5 3 (a;)| < eg (1 + | ln^|) 2 ||H||. 
Therefore this term is bounded like 

vWSzfm-L-M) <4 M (i + |in M |) 4 ||F|| 2 



for some constant Cg. 
We complete the proof adding these three inequalities. 



□ 



To pursue the analysis, we begin by rewriting the general form of the integral 
equation (|3.9p . showing that the various singularities of the equation can be 
rccombined under a more convenient form. This intermediate result is essential 
to obtain all following results. Indeed the integral equation for U (|3.9[) choosing 
G = writes 

(a(x) + ifi)U(x) = a V 9 A^(x) + b V x B fl (x) + [ V 9 x V 9 z k tl {x,z)U{z)dz. 

Jo 

Since by construction aoT> x A^(0) + boD x B^{Q) = Q(U) one also has 

(a(x) + in)U(x) = a Q (V 9 x A^x) - 2^(0)) + b (V 9 B^(x) - V 9 B^0)) 

+Q(U)+ f V 9 x V 9 z k^{x,z)U{z)dz. 
Jo 

But one also has due to the integral equation for V (|3.12j) choosing G = H 
V(0) = a = anA^O) + bnB^O) - / T> e z k»(0, z)U{z)dz. 



Basic manipulations yield 

a-o = an - 

Jo 



(V 9 z k^{0, z) - V 9 z k^ l (0 7 0)) U{z)dz - id \ U(z)dz 



H 
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because 2?^fc M (0,0) = i6. Since the function T> 9 z k^ is continuous, there exists a 
constant C% independent of \i such that 

\V e z k^{x,z)-V e z k^{x,x)\ < Cf(z-x) < Cfz for < x < z < H. 

Therefore the integral 

H rH 

\V e z k'"(0,z)-V 9 z k fi (0,0)\\U(z)\dz<Ci / z\U(z)\dz 

Jo 

is bounded uniformly with respect to /j, thanks to the bound given in (|5.20l) . We 
summarize this as 

f H 

ao = a — i9 / U(z)dz (5.31) 
Jo 

where \a\ < Cf H-ffH is bounded uniformly with respect to \i. Similarly 

b = b H - d x V e z k fl (0,z)U(z)dz (5.32) 
Jo 

and since the function d x T> z k^ is continuous and d x T> z (0, 0) = £5(0) 

b = b-i / 5(0)U e ^(z)dz (5.33) 
Jo 

where b is also bounded uniformly with respect to fi: \b\ < C 9 \\H\\. The integral 
equation then gives 

(a(x) + in)U{x) = a (V^A^x) - T> X A^Q)) + b (v£B M (x) - V e x B^Q)) 

+Q(U)-/ Q(x,z)U(z)dz + f V x V e z V L (x,z)U{z)dz 
Jo Jo 

where the new kernel is 

Q(x, z) = (V x A^x) - 2?«A„(0)) i6 + (x£fl M (x) - 2>£fl„(0)) iS{z) 

= T> x A»(x)V e z B»(Q) - V e x B»{x)V e z A»(Q) = V x V e z k»{x, 0) 

after evident simplifications. It is convenient to introduce two bounded functions 

m ^ = — « — — s — — and n ,M = — — — - — — so that (|3.Q|) is rewritten 

as 

(a(x) + %n)U{x) = ~am e ^(x)x + bn e ^(x)x + Q(U) - / V 9 x V e z k^{x, 0)U(z)dz 

J X 

(5.34) 

+ / (V e x V e z k^(x,z)-V e x V e z k fl {x,0))U(z)dz, Vxe[-L,oo[. 
Jo 

A first property which shows that (|5.34[) is less singular that its initial form 
(|3.9p is the following lemma which uses the pointwise estimate (|5.20[) on U (so 
an important restriction is nevertheless that x > 0). 
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Lemma 5.2. The first component U of any element U € X ,M satisfies 



(a(x) + ifi)U(x) = p°^(x)x + Q(U) 
where 



V e x V° z k fl (x,Q)U(z)dz (5.35) 



llp^lU-fo^^C^IIFH, V/i£ [0,1]. 



(5.36) 



Proof. Let us focus on the second integral in (|5.34[) . Continuity properties 
with respect to the second variable z imply that there exists a constant C 7 
independent of fi such that 



\V e x D e X{x,z)-V e x D e z k»{xM < C 9 7 z. 



(5.37) 



< C e 7 / z\U{z)\dz< C e 7 C e \\H\\x 



So, for x > 0, 

(V^V^ix, z) - V e x V e z k^{x, 0)) £/(z)dz 
using estimate (|5.20[) . Set 

p e '"(x) = ~am e ^{x) + bn B ^{x) + - f {V e x V e z k^{x, z) - V e x V e z k"(x, 0)) U{z)dz 

x Jo 

(5.38) 

which satisfies by construction (|5.36[) . □ 

As a consequence one has 

Proposition 5.12. For all 1 < p < oo, there exists a constant C® independent 
of /i and which depends continuously on 8 such that 



U 



Q(U) 



a(-) + iyU 

Proof. From lemma 15.21 one has that 



L"(0,H) 



<C e p \\H\\. 



(5.39) 



U(x) 



Q(U) 



V s x V e z k»(x,0) f H 



a(x) + i/i a(x) + z/i 
which turns into 

U(x) 



a(x) + ifi 



U {z)dz, 



a{x) + ifi J a(x) + ifi 



U(z)- 



Q(U) 

a(z) + ifi 



dz 



= ——-p e ^(x)-Q(V) 
a[x) + ifi 

By virtue of ()H2[) we notice that 



V e x Vl^(x,0) f H 1 



a(x) + i/i Jj. a(z) + ifi 



-dz 



(5.40) 



1 



a(z)+z'/i 



f?2 



< 



H 



\a{zj\ dZ -^ 0g{HM - 
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Since all powers of the function x i-> In \x\ are integrable, the right-hand side 
(|5.40[) is naturally bounded in any L p , 1 < p < oo. Therefore the function 
Z(x) = U (x) — * s solution of an integral equation with a bounded kernel 

and a right hand side in LP . The form of this integral equation is 

Z{x)+K 8 ' fl {x) / Z{z)dz = b e '»{x) 

J X 

with K e ' tJ, (x) < Cf independently of One also uses H^'^H^q ff -> < 

Cp 1 1 77 1 1 for < /i < 1: the key estimate is (|5.36j) which explains why the result 
is restricted to x > . Since this is a standard non-singular integral equation, 
see [55], the claim is proved. □ 

The previous result (|5.39j) shows that some singularities of the integral equa- 
tion can be recombined in a less singular formulation, so that the dominant part 
of U is a (.) +i/i ■ An important restriction of this technique, for the moment, is 
that it needs the a priori estimate (|5.20[) on U. This explains why inequality 
(|5.39[) is restricted to x > 0. By inspection of the structure of the algebra, 
it appears that one has the same kind of inequalities on the entire interval by 
replacing U directly by the function a ^ +ifJ in the integrals. A preliminary and 
fundamental result in this direction concerns the function 

D e ^{x) = x .\ \ ' ——^ -dz+ x \ V —r\ — ~ dz 

a{x)+tfi J a{z) + ifi J a(x) + i\i a(z) + i\i 

which is nothing than the integral part of (|5.34j) where U is replaced by the 
function , } , . . 

Proposition 5.13. Let 1 < p < oo. One has \\D 8,IJ '\\ LP ^_ L h ^ < C 8 where the 
constant depends continuously on 6 and does not depend on p. 

Proof. Two cases occur. 

• Assume < x < H . The analysis is similar to the one of proposition 
15.121 One has the same kind of rearrangement (|5.34[) , that is 

Ztplpx^ I* 1 ^ 
a[x) + ifj, J x a(z) + in 

V 8 V 8 k> i {x,z)-V x V e z kv(x,0) 1 

a(x) + i\i a(z) + ifx 

The first term is bounded like C el -^1 which is in all L p , p < oo. The 
second term is immediately bounded using (|5 .37|) : indeed 



Jo 



V 8 x V 8 k»(x, z) - V 8 V 8 z k»(x, 0) 1 ^ 
a(x) + i\i ol(z) + ifx 

< c 8 - 1 r - z dz < c 8 \. 

- 7 yju{xf + pi 2 Jo y/a{ Z ) 2 +fl 2 - V 
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• Assume — L < x < 0. The decomposition is slightly different and uses 
some cancellations permitted by the symmetry properties of the kernels. 
One has 



D e ^(x) = - 



H 



1 



V x V 9 k^(x,0) 
a(x) + ifi J_ x a(z) + ifi 



-dz 



x T,e 



V x V°k»(x,z) 



1 



o u(x) + ifj, a(z) + in 



dz — 



V x V 9 z V L (x,Q) 



1 



a(x) + ifj Jo a{z) + t/i 



dz, 



which emphasizes the importance of some symmetry properties of the 
kernels. Indeed 



o a(w) + ifi 



o u(z) + ifj 
dw 



dz 



1 



-dw 



1 



o a(— wj + iix 
1 1 



dw. 



Notice that 



1 



-a{— w) — ifi a(w) + ifj 

a(w) + a(—w) + 2i/j 
— a(— w) — ifj a(w)+i/j (a(w) + ifi)(—a(—w) — ifi) 

So, since a(0) = 0, 

1 1 



< 



2 W a \\w^(-LM) W 2 + 2/1 



-a{— w) — ifj a(w) + ifj 
since a £ W 2 '°°(—L, H). One can bound 

1 



r 2 w 2 + fi A 



-dw — 



o — ct\— w) — ifi 



1 



< 



< 



\W 2 -°°(-L,H) 



o a(w) + ifi 
2/t 



dw 



r 2 z 2 + fi l 



dz 



H 



2fi 



r 2 z 2 + fi z 



T dz < 



\W 2 -°°{~L,H) 



H 



7T 

r 



As a consequence D 9 '^ can be expressed as 

V 9 V 9 k^(x,0) f H 1 



D 9 '"(x) 



a(x) + ifi J _ x a(z) + ifj 
x ViV<> z W{x,z)-ViVlW{x,Q) 1 



-dz 



a{x) + ifi 



a(z) + ifi 



dz + R(x) 



with H-^Hoo (— L, H) < Cf . The two integrals have the same structure as 
for the first case, in particular the interval of integration is [— x, H] with 
< —x. So the same result holds. 
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□ 

Proposition 5.14. For all 1 < p < oo, there exists a constant C 9 independent 
of /i such that 

Q(U) 



u 



a(-) + ifj, 

Proof. We start from (|5.34|) written as 

Q(U) 



Lp(-L,H) 



<C;\\H\\. (5.41) 



U{x) 



a(x) + ifi a(x) + ifi 



p°^(x) 



U X U K [X,U) u , z)dz _ 

a[x) + ifi j 



V e x V e z W{x,z) 
a{x) + ifi 



U(z)dz. 



Here p e ^{x) = hm 9 ^{x) + bn e ^(x), so that Wf^W^i-LM) ^ C<> \\ H \\ over the 
whole interval (-L, H). Notice that p 9 ^ is the first part oip '^ defined in (|5.38l) . 
Setting u{x) = U{x) - a f}^ one gets 

U(X) — / t— r ; U(Z)dZ 



o a{x) + ifi 



* ^(x)-Q(U)D 8 >»(x)- 
a(x) + ijj, j 



H V s x V e z k»(x,0) 



u{z)dz. 



a(x) + ifi 

The left-hand side is an non singular integral operator of the second kind with 
and with a bounded kernel thanks to the fundamental property (|4.4[) . The right- 
hand side is bounded in LP with a continuous dependence with respect to 
see Lemma [521 estimation (|5.29|) and estimation (|5.39j) . 

□ 



5.4.4 The second basis function 

We apply the above material to the second basis function for which Q(U 2 ' M ) = 1. 
The inequality (|5.41[) writes 



U, 



1 



a(-) + ifj, 
for 1 < p < oo. 



Lp(-L,H) 



<C e p (\v 2 e ^(H)\ + \w^(H)\), (5.42) 



Proposition 5.15. Assume the transversality condition i5.19\) . There exists a 
constant D e independent of jj and continuous with respect to 9 such that 



Vf'^H) + W^{H) < C 



(5.43) 
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Proof. Indeed, regarding relation (l5~Ti)) . (15~T5l the pair (v, w) = {V^{H), W° 2 -"(H)) 
is solution of the linear system 



-vW^{H) + wV°^{H) = 1, 
vW^{H) - wVi'"(H) = 0. 



(5.44) 



The determinant of this linear system is equal to the value of the function 
—a(9,fi). So the transversality condition establishes that 



det 



-Wf*(H) V^(H) 
W^{H) -V^{H) 

Therefore the solution of the linear system 



is bounded uniformly with respect to ji. 



(5.45) 



□ 



Theorem 5.1. Assume the same transversality condition. The second basis 
function satisfies the following estimates for some C e v and C e which are contin- 
uous with respect to 8 



U. 



1 



«(•) + ifi 



Lp{-L,H) 



<C°, l<p<oo, 



-L,0)U(0,H] 



(5.46) 



(5.47) 



Proof. The first estimate is a straightforward consequence of (|5.42[) , (|5.43l) 
The use of the integral representations (|3. 12113. 13")) shows that, 



O./i 



-L,0)U(0,H] 



-L,0)U(0,H] 



(5.48) 



for some C e . Then the second equation of 
bound for U. 



O41 



u 2 



L~ c [-L,0)U(0,ff] 

The bound on the derivatives follows from 



shows that one has the same 

< C°. (5.49) 

□ 

Remark 7. Let us set W = —L. From Hk^M one gets that \\H'\\ is bounded 
uniformly also, therefore 115. 20\) can be generalized for x < ( resp. H' ) instead 
of x > (resp. H). In summary one has for a constant K e that can be further 
specified 



U^(x) 



< 



9 9 1 9 : 
r z x + // 

We now pass to the limit /i — > + . 



x e (-L,H). 



35 



Proposition 5.16. Assume the transversality condition. The second basis func- 
tion admits a limit in the sense of distribution for fj, = 0^ as follows: 



U a '"->U 9 =[py—±—8 D + u 2 ,v 2 ,w 2 



where u 2 '^ , v 2 '^ , G L 2 {— L, oo) and Sd is the Dirac mass at the origin. 
Remark 8. The limits are solutions of \2. 3]) in the sense of distribution. 



they will be called the singular solutions. 

Proof. We consider the case fj, I 0. Some parts of the proof are already evident, 
essentially for quantities which are regular enough (V^ 8 '^ and W 2 ,/< ) or for regions 
where all functions are regular (typically x > 0). Therefore the whole point is 
to pass to the limit in the singular part of the solution U 2 ■ We will make wide 



use of the equivalence between the integral formulation of proposition 13.11 and 
the differential formulation (|2.3j) . 

• Passing to the weak limit: By continuity of the first basis function with 
respect to fi, one can pass to the limit concerning (V 2 '^(H), W 2 '^(H)). One gets 
that (v,w) = (V 2 '° (H), W 2 '° (H)) is the unique solution of the linear system 

-vW?(H)+wV 1 e (H) = l, 
vW$(H) - wV£(H) = 0, 

where the coefficients are defined in terms of the first basis function for /i = 0. 
By continuity away from the singularity at x = 0, one has that — > U2 in 
L°°{e, H) for all e > 0. Using ([5~3o]) it is clear that U '^ - a( .) +lfJi is bounded in 
L 2 (—L, H) uniformly with respect to ^. Therefore there exists a limit function 
denoted as u 2 such that for a subsequence 

U 2 r— . > W eak U 2 ' m L (-L, H). 

Oi(-) + Ifl 

Moreover the first derivative of U 2 is bounded in L 2 (—L,—e) by virtue of 
([5~17|) . Therefore 

U 2 ' P ^strong — TT + U 2 ' in L {—L, —£) 

a(-) 

at least for a subsequence. Considering the integral relations (|3. 12113. 13"|) . these 
subsequences are such that 

V 2 e ^(x)->v 2 O+ (x), (5.51) 

and 

W°>»{x) ^ w°'° (x), (5.52) 
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with the convergence uniform in compact sets of (-L, H) \ {0} and where, 



9,0+ 



a H A (x) + b H B (x) + f V x (k\x, z) - k°(x, 0)) ( + u e '+(z)) + 
V x k°(x,0)u e -+(z)+d(x), 



with 



H 



v(x):= V e x k°(x,0) -^dz, forxX), 



lna(-e)/3(-e) 



v(x) :=V e x k°(x,0) 

,a(-e) 



1 

a(z) 
1 

x Ot{z) 



dz + lna(e) /3(e) - 



a(e) 



ln(2) /3'(z)dz 



Vlk"{x,Q)——dz for,x<0, 
a(z) 



w e 2 fi+ := a H A' (x) + b H B' (x) + j' d x V e x (k°(x, z) ~ k°(x, 0)) ( J-y + u e > + (z)^j 



+ I d x V x k°(x,0)u 9 ' + (z)+w(x), 
'h 



with 



X 



1 



w(x):= d x V x k u {x,0)— r ^dz, forx > 0, 



H 

6 1.0/ 



w(x) := d x T> x k u (x,0) 

r a(-e) 



a(z) 
- e 1 



x a ( z ) 

H 



dz + \na(e) j3(e)- 



\na{-e)/3(-e) + / \n{z) (3' (z) dz 



1 



d x Vlk»(x, 0) — - dz for, x < 0, 
31 v ; a(z) 



where < e is so small that a(x) is invertible for x <G [— e,e] and f3(z) := 
l/a'(a~ 1 (z)) with a -1 the inverse function of a. The limits in (|5.51[) . (|5.52[) 
also hold in the strong topology of L 2 (~L, H). 

These weak or strong limits are naturally weak solution of the initial system 
(|2.3I) : denoting for simplicity 

/ \ / e.o+ e,o+ e,o+\ 

(U 2 ,V 2 ,W 2 ) = [U 2 ' ,V 2 ' ,W 2 ' ), 



these functions arc solutions of 



w 2 ipidx + i6 P.V. 



1 



a(x) 



On 



+ u 2 ) ipidx - ^7^yVi(0) + / v 2 <p\dx = 0, 



i6 J w 2 ip 2 dx — J (au 2 + l)ip 2 dx — i J S(x)v 2 ip 2 dx = 0, 

wttidx + i P.V. J S + u 2 ^j V ,dx - ^^3(0) 

— J a{x)v 2 (pzdx = 0, 



(5.53) 



37 



for any sufficiently smooth test functions with compact support, for example 
{(pi, v?2, V3) € L, H). To pass to the limit we have used that in distribution 



P.V- 



Thc signs of -^^(0) and 
^, AT 1^3(0) are compatibltH with the fact the limit is for positive ji. The 



sense, lim 

(O)tt 

m 



i compatible 
principal value is defined as 

1 



P.V. 



f(x) dx := lim 
a(x) ej.0 



-L Ol(x) 



Cf(x) 



1 



(p(x) dx, 



where a(p(^fej) = ±e. 

• Uniqueness of the weak limit: If there is another triplet (u2,V2,W2) 
solution of the same weak formulation (|5.53j) . then the difference 



(u 2 ,V 2 ,W 2 ) = (u 2 - U 2 ,V 2 - V 2 ,W 2 - W 2 ) 



satisfies 



wiipidx + i9 J u^fidx + J v 2 (p' 1 dx = §, 
9 j w^ip^x — J aui^dx — i J 5(x)v 2 (p 2 dx = 0, 
w^ip^dx + i / Su^ipzdx — / a(x)v 2 (fi3dx = 0, 



(5.54) 



By construction (wj, v^iWi) = (0,0,0) for x > 0. For x < 0, we deduce from 
(|5.54[) that (u2, v~ 2 , wi) is a solution of the X-mode equations. Therefore these 
functions can be expressed as a linear combination of the first and second basis 
functions for x < 0. Since e L 2 {— i,0) is non singular, only the first basis 
function is involved that is 

From (|5.54j) we get for example 



— L 



w^ip^dx + i / Su^ipzdx — / a^v^tps dx = 



-L 



-£ 



where ip^(-L) = and </93(0) is arbitrary. We integrate by parts 



(— W2 + iSui — av 2 ) (f 3 dx + £2(0)^3(0) = 0. 



Since iui, v 2 , w 2 ) is a non singular solution of the X-mode equations, one has 
that — w^' + iSu^ — av 2 = 0. Finally 

£E(ov 3 (o) = 0. 



J If one takes the limit fi 0, the signs of these terms are changed 
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Since we can take ^3(0) 7^ 0, it follows that = ^(0) = XWf (0). Considering 
the normalization (|5.6[) one gets that A = 0. Therefore (wj, v~2,W2) = (0,0,0). 
It means that the weak limit is unique: all the sequence tends to the same weak 
limit. 

• Regularity: By Theorem 1 5. II the limit belongs to H 1 {[—L, — e] U [e, oo)) 3 . 

□ 

6 The limit spaces X 6 ^ 

We can now define the limit spaces in which the limit basis functions live. 

6.1 The space X e > + 

Passing to the limit \i — > + , the limit space X e,+ is 

Xf ■+ = Span {Uf, U*<+} C U\ oc ((-L, 00) \ {0}) . (6.1) 

6.2 The space X 9 ~ 

It is of course possible do all the analysis with negative fi < and to study the 
limit fi — > 0~. The first basis function is exactly the same. The second basis 
function is chosen exponentially decreasing at infinity and such that 

i/jtU^ = 1 n < 0. 
The generalization of the preliminary result (|5.46[) is straightforward 

< C e pl -1 < /x < 0. (6.2) 

LP(-L,H) 

Passing to the limit /i — >• 0~, it defines the limit space X e '~ 

= Span {U?, U|'- } C H\ oc ((-L, cx>) \ {0}) . (6.3) 

We observe of course that the first basis function belongs to X^ ,+ n Xf ,_ . Since 
the limit equation is the same, and the normalization at x = H is also the same 
we readily observe that the second basis functions are identical for < x 

U^ + {x) =U%'~(x) 0<x. (6.4) 

The main point is to determine the difference between the limit of the two 
singular functions for x < 0. 

Proposition 6.1. One has 

U 9 2 + {x) ~ U e 2 '-(x) = ^Ul{x) x<0. (6.5) 



U^'* 1 TV — ~ 

a(-) + t/i 
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Proof. We notice that the Wronskian relations (|5 . 1 5[) are the same at the limit 
[i = 0^. By subtraction 

V?(x) (w 6 2 + {x) <») - W°(x) (v 2 e ^(x) - V 2 e '~(x)) = 0. 

It show that the difference is proportional to the first basis function 

U 2 + (x) - U 2 '~(x) = 7U1 (as) x<0. (6.6) 

It remains to determine 7. We already now that the limit fj, — > + can be 
characterized by (|5.53|) . The third equation writes 

J w+(p' 3 dx + i P.V. J 6 + u t^j ¥3dx - ^^3(0) - J a{x)v^tp 3 dx = 

where (u 2 ,v 2 ,w 2 ) refers to the non singular part of the limit /i — > + . The 
equivalent equation for the non singular part (u 2 ,v 2 , w 2 ) of the limit fj, — > 0~ 
is 

w 2 ip' 3 dx + i P.V. J 5 (^j^j + "2 J fsdx + 13 (0) - J a(x)v 2 ip 3 dx = 0. 
By subtraction, one gets 

iw 2 — w 2 )ip' 3 dx + i J S(u^ — u 2 )ip 3 dx 

-^^y-^3(0) - J a{x)(vt - v~)ip 3 dx = 0. 
Due to (|6.4[) these differences vanishes for x > 0. We get 



(w 2 — w 2 )(p' 3 dx + i I S(u 2 — u 2 )tp 3 dx 

J-L 

-^—^3(0) - / a(x)(v+ - v 2 )cp s dx = 

"(0) J-L 

where ip 3 is a smooth test function that vanishes at —L. Integration by part 
yields 

r 

{-( w t ~ W 2 Y + i5[ui - «2 ) - a ( x )( v 2 ~ v 2)) <P3 dx 

2<y(o)7r_ ^ , u + 



a'(0) 

Due to (16.61) one has that 



p 3 (0) + (w+ - w 2 )(0)ip 3 (0). 



~(w+ -w 2 )' + i5{u+ -u 2 )- a(x)(v£ -v 2 ) = x < 0. 

Since tp 3 (0) is arbitrary, it means that w 2 (0) — w 2 (0) = 2 ^,?q\ ■ We obtain 

7< (0) = that is t S(0h = ^gf = 0. Therefore 7 = The claim 

is proved. □ 
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7 Proof of the main theorem 



All the information about the first and second basis functions is now used to con- 
struct the solution of the system (|2.3[) with the boundary condition (|1.4p . The 
function g depends only of the vertical variable y. Under convenient condition 
g admits the Fourier representation 




We first consider a small but non zero regularization parameter fi > 0. For the 
sake of simplicity we will assume that the transversality condition is satisfied 
for all 9 in the support of <? 

|ct(0)| > c> V0 e supp (g) . (H6) 

It is just a convenient uniform version of the point-wise transversality condition 
(|5.19p . Additional comments are in section [7731 



7.1 One Fourier mode 

For one Fourier mode, one needs to consider the solution of (|2.3p with boundary 
condition 

W^-L) +iXV tl (-L) =g. 

Since we add of course that the solution must decrease (exponentially) at x ~ oo 
to guarantee that no energy comes from infinity, the solution is proportional to 
the second basis function. That is there is a coefficient j 9 ^ such that = 
ry6,t*\j6,n^ The coefficient satisfies the equation 

7 ^ (W*'»(-L) + iW^(-L)) = g{6) 

that is 

~m = m 

from which it is clear that we must study the coefficient / function 



r 



fl,M = W^{-L) + iXV^i-L). (7.2) 



Proposition 7.1. Assume IIH6\) . For every compact set S C 1, there exists 
e > 0, t + and t_ > such that t~ < \t 9 ^\ < t+ for < (i < e and 9 e S. 

Proof. The upper bound is a direct consequence of (|5.48p . To prove the lower 
bound, a useful result is the formula which comes from (|5.10p 



Im 



{W^(-L)V 2 e ' v (-L)) >»\ l U^ix) 



2 

dx 
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Combining with (|5.27[) and Q (U 6 ^) = 1 (by construction), it yields 
Im (w^{-L)V^{-L)) > r_ > 0. 



Plugging the definition of r ,fl inside this inequality, one gets 

Im (i- fl -"Vf '"(-£)) > r_ + A \V e '^{-L)\ 2 > r_ > 0. 

Therefore V^{-L) x |r(0,/i)| > r_. The i°° bounds (|5^5|) shows that there 
exists C > such that C |r(0, fi)\ > r_. □ 

By (EHD, (O, 

r e <+ := ^ 2 M+ (-L) +*A^' 0+ (-L) = lim r^. (7.3) 

Proposition 7.2. For every compact set S C I, there exists e > 0, t + and 
r_ > such that r~ < |r e ' + | < t+ /or < \i < e and 9 G 5. 



Proof. This is immediate from Proposition 17. II and (|T.3|) . □ 
7.2 Fourier representation of the solution 

The solution of (|2.2I) with the boundary condition (|1.4[) is given by the inverse 
Fourier formula 



\^y) = ^i ^U e ^(x)e iey de (7.4) 



where we assume that g <E L 2 (R) and that g has compact support. Passing to 
the limit in (|7.4[) one gets 



i /■?(*) 



V. 



2 



(7.5) 



This formula has to be understood in the sense of distribution. Since by Theorem 
EZD||u 9 ' + H < C|) |K' + j| ^ C 2>|| u ' e ' + || < ^2 withCf a continuous function of 9, 
and considering that r ,A1 converges to r e,+ there is sufficient regularity to pass 
to the limit. The value of the heating is 

a = ( to „/ Kfe »)!■** - i / r jjj^j,* (") 

To our knowledge this is the first time that such a formula is written where all 
terms are explicitly given. A similar but much less precise formula can be found 
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in [TT] derived by means of analogies, see also [55]. All these integrals are conver- 
gent provided g decays sufficiently fast at infinity. Provided the transversality 
condition is satisfied, and it is always the case in a neighborhood of 6 = under 
hypothesis (|H5|) . the heating is generically positive. This is of course related to 
the presence of the strong Dirac singularity in the solution (|7.5[) . 

Remark 9. An essential consequence of this analysis is the physical heating 
Q which is related to the singularity P-V.-^^ ± a " ^ <5p of the mathematical 
solution. The singularity is not an artifact of the model. It is on the contrary 
a direct way to measure the amount of heating provided to the plasma by the 
electromagnetic wave. 

Remark 10. Observe that the singular solutions U 2 ' are the unique solutions 
of the following initial value problem. 

Find a triplet (u!^, u, 2' ± ) G L 2 {— L, oo) 3 which satisfies the constraints 
= V£'°(H), w^iH) = Wi'°(H), and 

W2 __ W2 . +l e U2 = _^p.y.__±__ fo) 

■a o,± e,± -x? \ e,± -, 
it)w 2 — au 2 — io{x)v 2 =1, 

d s.± . ■ e,± / n e,± • D T/ 5 ( a; ) . ^(^x 
— — w 2 +iu 2 —a[x)v 2 =-iP.V.——± .,,., o D . 
dx 1 * * a(x) a'(0) 

We prove that this problem has an unique solution by the argument given to 
prove the uniqueness of the weak limits.. For this purpose observe that 

(^y +t4' ± l t£ ± ,«£ ± ) (*) = (u!°,Vr,Wl°) (x) forx> 0. 

We observe the similarity with the standard limiting absorption principle in 
scattering theory. In scattering theory the solutions obtained by the limiting 
absorption principle are characterized as the unique solutions that satisfy the 
radiation condition, i.e., they are uniquely determined by the behavior at infinity. 
Here, the singular solutions are uniquely determined by their behavior at +oo 
and by their singular part P.V.-^^r ± Sd Note that it is natural that we 
have to specify the singularity at x — because our equations are degenerate at 
x = 0. We think this principle could be used for practical computations. It is 
however a little more subtle since a boundary condition at finite distance x = —L 
must be prescribed. That is the singular part is itself dependant on the boundary 
condition where the energy comes in the system. Mathematically it corresponds 
to the coefficient t 9 ' + in the representation formula |7.5[ ). 

7.3 What happens if the transversality condition is not 
satisfied 

An interesting question is to determine what happens if the transversality condi- 
tion is not satisfied. We begin with some simple remarks. Firstly the point-wise 
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transversality condition (|5.19|) or the uniform one (|H6|) greatly simplify the anal- 
ysis. They are satisfied at least for 9 close to zero provided the transition zone is 
small enough (|H5|) . Secondly some technical intermediate results may be wrong 
if these conditions are not satisfied: for example it is not clear whether a is still 
regular at points 9 such that a (6) = 0. 

Our purpose is not to answer to the questions raised by this possibilities, 
but only to give understanding of the physical situation hidden behind and to 
explain what is the limit value of the heating (|7.6p . The analysis is as follows. 

Physical picture: If a(9) = then the first basis function is propor- 
tional to U3 (at least for x > H). It means that it is also exponentially decaying 
at infinity. That is the excitation provided by the boundary condition "captures" 
this non singular first basis function which is also the physical one. We remark 
that no heating is provided by Ui because Q(Ui) = by definition. Therefore 
a shortcut is: if a(9) = 0, then physical heating vanishes. 

Mathematical picture: On the other hand it is also clear that the function 
t°- + is upper bounded by virtue of the analysis provided in section 17. 1 1 A 
consequence of the transversality condition is the fact that is uniformly 



In 



bounded, see (|5.48|) . So if a (6) = 0, it is possible that lim^o 

this case \t{9)\ = 00 which yields in turn once again that the associated heating 
vanishes in (|7.6p . 

In summary it is possible to conjecture that (|7.5[) is still valid even if the 
transversality condition is wrong: in this case the heating associated to the 
Fourier mode vanishes. The limit of (|7.4[) may a priori be more singular. More 
research is nevertheless needed to provide a rigorous basis to this analysis. 

A Appendix 

We provide additional material to obtain a self contained paper. 
A.l Approximation of Airy functions 

Suppose that A and B are the two fundamental solutions of the equation — u" — 
au = satisfying the normalization conditions 



A(0) = 1, A'(0)=0, 
B(0) = 0, B'(0) = I, 



(A.1) 



such that the corresponding Wronskian is equal to 1. Let A^ be an approxima- 
tion of A in the following sense 



-A" - aAf, = f Afi , with f Ait := ifxA^, 



11 

^(0) = 1, ^(0) = 
From the variation of constants one gets 



A^x) = A(x) (c A + J f A ,(t)B(t) dt\ + B(x) L B - J f A ,(t)A(t) 



(A.2) 



dt 
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The initial values (|A.2[) yield A^(x) = A(x)+ /a m {t)k{x, t) dt, where k(x, y) = 
A(x)B(y) — A(y)B(x). So A^ satisfies a classical Volterra integral equation 

A^x) -in [ Af,{t)k{x, t) dt = A(x). (A.3) 
Jo 

Define the series of integral kernels 

f K 1 (x,y) = k(x,y), 

\ K n+1 (x,y) = Jo k(x,x n )K n (x n ,y) dx n . 
The solution of the integral equation (I A. 3 1 is 

A„(x) = A(x) + 

For n > 1 

*•.„<*.»)-/ n *m.*> *«*<*.») to, 

J0<Il<~<I„<I i^.^„ 



l<i<n 

and 

^n(- T ) = J( xi ,...,a;„)e{0<a;i<'-'<x n <x} ril<i<n 
— J*^ ^n— l(*^n) dx ni 

= ^-f since ii(ir) = x. 
So the iterated kernels satisfy > 

|JT Tl+1 (x > »)|<(2||A|| 00 ||B||J"+ 1 x7n!. 

On the compact interval ]0,L+[ the sum and integral symbols can be inverted, 
which gives with a shift of the index n 



A,(x)=A(x) + Y / [ / 

n=l VO 



x,y)A(y) dy ) (i/x) 



Assuming is bounded positive for the simplicity of notations, A^ is indeed 
bounded independently for < fi < 1 

with C = 2^+11^11^1151100. From (|X3|) it yields 

Similarly if approximates B in the following sense 

—Bp - aB v = ifJ-Bft, 
fl M (0) = 0, fl;(0) = l, 
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then one has the inequality < H-BH^ e.^ Ca together with 

\B^x)-B{x)\< i iC^\\B\\ oa e^\ 

Both Ajj_(x) and B^x) are C°° functions with respect to \i and x. Since any H 1 
function / is 1/2 Holder thanks to the inequality \ f(x) — f(y)\ < \\P\\l 2 \ x ~ 2/| , 
then Afj, (x) , B^ (x) as well as all their derivatives with respect to \x and x are also 
1/2 Holder, with constants bounded independently of fi as far as fi is bounded. 



B Derivation of the cold plasma dielectric ten- 
sor 



Our aim is at the derivation of the cold plasma dielectric tensor used in this 
work. Much more material can be found in classical physical textbooks [161 [7] . 
Let us consider Maxwell's equations with a linear current 

-49 f E + VAB 



<9 t B + VAE = 0, 

m e dtU e = — e (E + u e A Bq) — m e vu e 



-eN e u e 



(B.l) 



Here the linearized equation for the electrons corresponds to moving electrons 
(velocity u e ) with a given electronic density N e which is given but may be non 
constant in space. One implicitly assumes that an ion bath, which is the reason 
of the friction between the electrons and the ions with collision frequency v. 
Shifting to the frequency domain (that dt = ~iuj) yields 

jziuF, + V A B = -/j. eN e u e , 
-itoB + V A E = 0, 

—im e Luu e = — e (E + u e A Bo) — m e v\x e . 

One can compute the velocity using the third equation 

g 

um e + Lo c i\i e A bo = iE 



(B.2) 



(B.3) 



where the cyclotron frequency is uj c — ^— - ■ , bo — JSoT ^ 

the normalized mag- 



^,b 

netic field and uj = oj + iv is the equivalent a priori complex pulsation. This is 
a linear equation. Assuming that bo = (0,0, 1) one gets 



u, 



V 



UJ* —to* 





E. 



—i 



UlrCO 



WW 



\ 



o 












(B.4) 



It is then easy to eliminate u e from the first equation of the system (|B.2[) and 
to obtain the Maxwell's equation (|1.1[) with the dielectric tensor 



(B.5) 
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The dielectric tensor is slightly more general than (|1.2p in the sense that uj ^ uj 
a priori. However physical values in fusion plasmas are such that the collision 
frequency is much smaller than the frequency (y « uj) which means that some 
simplifications can be done in the dielectric tensor. We refer for example to |16j . 
page 197. In this case one recovers exactly the limit tensor (|1.2[) . This derivation 
is the explanation of the non zero extra-diagonal terms in the dielectric tensor. 
Notice that we consider in this work uj 7^ u> c , that is the frequency is away of 

the cyclotron frequency. Furthermore, when v ^ the diagonal components of 

~ 2 

the dielectric tensor, 1 — ^(u 1 —^ 1 ) nave an imaginary part. This shows that our 
rcgularization parameter \x can be interpreted as a small collision frequency. 

The loss of energy in domain SI can easily be computed in the time domain 
starting from (|B.1[) . One obtains 



boundary terms. 




Therefore 



vrripNp I u„ 



S2 



represents the total loss of energy of the electromagnetic field plus the electrons. 
This loss of energy is necessarily equal to what is gained by the ions. This 
expression can be evaluated either in the time domain or also in the frequency 
domain. Considering (|B.3[) another equivalent formula in the frequency domain 
is 

2 = - Re 



/ eiV e (E,u e ) 
Jn 



Finally one can eliminate the electron velocity in function of the electric field 
using (|B.4[) . Therefore a third formula is 



Q = -cj£ Im 



/ (E,eE) 
Jn 



where e is the dielectric tensor (|B.5[) . These formulas are the physical analoguous 
to the mathematical formula 
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